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Abstract. We propose a scheme for precise control of tunneling dynamics of dipolar 
bosons in shaken triple-well potentials. In the high-frequency regimes and under the 
resonance conditions, we have analytically and numerically demonstrated that we can 
transport a priori prescribed number of dipolar bosons along different pathways and 
different directions by adjusting the driving parameters. These results extend the 
previous many-body selective coherent destruction of tunneling (CDT) schemes for 
nondipolar bosons in double-well potentials[Phys. Rev. Lett. 103, 133002 (2009); 
Phys. Rev. A 86, 044102 (2012)], thus offering an efficient way to design the long- 
range coherent quantum transportation. 
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1. Introduction 

In recent years, quantum-degenerate dipolar gases have attracted a great deal of 
attention from both theoretical and experimental studies|T]-|3]. In addition to the 
short-range and isotropic (s-wave) contact interaction, which is usually at work in 
ultra-cold gases, long-range and anisotropic dipole-dipole interaction (DDI) also plays a 
signihcant role in dipolar quantum gases and gives rise to a rich variety of new physical 
properties. A minimal system for direct visualization of the nonlocal characters of DDI 
is dipolar Bose-Einstein condensate (BEC) in a triple-well potential, which is described 
by a three-site Bose-Hubbard model with neighbor interactions. This minimal system 
has been actively studied and known to display some novel features, such as mesoscopic 
quantum superpositions[5], interaction-induced coherence [6], role of anisotropy [7] and 
entanglement entropy [8]. In addition, dipolar Bose-Einstein condensates in triple-well 
potentials have been also investigated by ways of mean-held treatments|9]-|TI] and multi- 
conhguration time-dependent Hartree (MCTDH) method[I^. 

Recently, there has been a burst of interest in experimental realization of quantum 
control with periodic lattice shaking technique. Periodically shaking lattice leads to 
many important physics, such as coherent destruction of tunneling (CDT)[T3l [HI 
ITS] , photon-assisted tunneling [16], Mott-insulator-superhuid transition[T7], simulation 
of frustrated classical magnetism|18] and effective ferromagnetism [19], control of 
cotunneling and superexchange [20] . generation of synthetic gauge held[2T], realization of 
Haldane model [22], and so on. Among these intriguing aspects, CDT effect is a simple 
and powerful tool to control quantum tunneling dynamics [231123], based on which two 
schemes for selective CDT of strongly interacting bosons in a symmetric double-well 
potential have been established by modulating the self-interaction strength [25] or energy 
level unbalance [26] . In these two many-body selective CDT schemes, the modulation 
can be tuned in such a way that only an arbitrarily and a priori prescribed number 
of bosons are allowed to tunnel from one well to the other. The sensitivity of CDT 
to particle number has also been demonstrated in an earlier work[27], which enables 
the self-trapped state to be used as a quantum beam splitter. As three-site system is 
a paradigmatic model for longer array of lattice, extensive efforts have been paid to 
study the dynamics of periodically shaken triple-well potentials [2H]-[3^. More recently, 
directed tunneling of 1 or A^ — 1 dipolar atoms with A^ > 2 in a triple-well potential 
has been found numerically [33]. It naturally leads to another question: is it possible to 
precisely control an arbitrarily prescribed number of dipolar bosons allowed to tunnel 
by shaking triple-well potential? 

In this article, we suggest a method to control directed tunneling of a prescribed 
number of dipolar bosons in periodically high-frequency shaken triple-well potentials, 
which extends the previous many-body selective CDT schemes for nondipolar bosons in 
double-well potentials[25l[26]. Under certain conditions, we analytically and numerically 
demonstrate that the driven three-site model can be reduced to a double-well model, 
in which only tunneling between two adjacent wells is allowed. We also reveal that 
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direct tunneling of a priori prescribed number of dipolar bosons occurs along different 
pathways and different directions. Our results may provide an additional possibility for 
designing the long-range coherent quantum transportation. 

2. Model and high-frequency approximation 

The starting point of our analysis is provided by a driven three-site Bose-Hubbard 
Hamiltonian in the presence of long-range DDI, which describes the tunneling dynamics 
of dipolar bosons in a shaken triple-well potential. Under the single-band tight-binding 
approximation, by expanding the bosonic held operator as 'ipij) = where 

0j(r) is the Wannier state localized at the jth site and aj is the corresponding atom 
annihilation operator, the three-site Bose-Hubbard Hamiltonian is given by[5l |33] 

H = Htun + .Hint + Hex 

with 
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Htun ^ ^ T 

H-mt = - 1) + Ui(nin2 + ^ 2713 ) -F U2nin3 

^ i=i 

Hex = £(t)(alai - a^aa), (1) 

where hj = is the particle number operator at the the jth site, v = 

f (ir0*(r)[— V^/2 -|- Utrap(i')]0*_ti(r) is the hopping (tunneling) rate between two 
sites, Uq = g J (ir|0*(r)|^ -f- / |0j(r)p0j(r')pVdd(r — r')drdr' characterizes the on-site 
interactions, Ui = f |0j(r)p0j_|_i(r')pVdd(r — r')drdr' is the coupling constant induced 
by nearest-neighbor DDI, and U 2 = J |0i(r)p(;/)3(r')pVdd(r — r')(ir(ir' describes the next- 
nearest-neighbor interaction. In the above expressions, the natural unit h = m = 1 
with m being the atomic mass is used, g = dvra^ is the familiar short-range interaction 
constant with being the s-wave scattering length, Utrap(i‘) is the time-independent 
triple-well scalar potential, and Vdd(i' — r') is the DDI potential. Here the external 
driving held is applied in the form[Tni [33], £{t) = Sq + Si cos{ut) with Eq and £1 being 
magnitudes of static and ac helds, u the driving frequency, respectively. The time- 
dependent part Hex of the Hamiltonian ([T]) can be experimentally realized by shaking 
(accelerating with a constant and modulating force) the triple-well potential|16j. 

The Hamiltonian ([1]) conserves the total number of particles. We can expand the 
vector state |'0(^)) of the system in Fock space according to 

W()>= E c„ \nun2,n3), (2) 

ni,712,713 

where rii -|- n 2 -|- 773 = N. Substitution of the ansatz ([2]) into the Schrodinger equation 
= H\'ip(t)), one obtains the following coupled equations for the probability 
amplitudes Cni,n 2 ,n 3 
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iCNflfi — VNflfiCNflfi + 

^N—n,n,oCj\f—n,n,0 4“ ^nC*JV—n—l,n+l,0 

4“ ^n— iC'tV— n+l,n—1,0 4“ l^nCpf—n,n—l,l 

(n= 

*C'o,Af,0 = ^,iV,oC'o,Ar,0 4* Kjsf-iCi^N-1,0 

+ I'^N-iCo^N-1,1 

iCo ,n,N-n — Vo,n,N-n Co ,n,N—n 4“ KnCo^n+l,N—n—l 

4“ l^n—lCo^n—l,N—n+l 4“ VyiC-i yi—i N—n 
(n = l,2,...,iV-l) 
iCofi^N = Vo,o,AfC'o, 0 ,iV 4- KoCo,l,N-l 

. J 

(3) 

where we have set 

,n2,n3 = (ni,n2,n3j(IIint + Bex)jni,n2,n3} 

Kn = — V\J {n + 1){N — n),n = 0,1, 2..., 

Pn = — v\/n,n = 1,2..., N 

VN-n,n,o = + eiCos{ut){N - n),n = 0,1,2, ...,N 

= y[(iV-n)(iV-n-1)+n(n-1)] 

4- Ui{N — n)n + eo{N — n) 

Vo,n,N-n = Fl - eicos{ut){N - n),n = 0,l,2,...,iV 

Fl = y[(iV-n)(iV-n-1)+n(n-1)] 

+ Ui{N — n)n — eo{N — n) 

VN-n,n-i,i = F^ + eiCos{u:t){N - n - l),n = 1,2,...,N 
Fl =^[{N-n){N-n-l) + {n-l){n- 2)] 

+ [/i[(iV-n)(n-l) + (n-l)] 

+ U 2 {N-n)Feo{N-n-1) 

Vi^n-i,N-n = F^-Si cos{ujt){N - u - 1),n = 1,2,..., N 

F^ =^[{N-n){N-n-l) + {n-l){n- 2)] 

+ [/i[(iV-n)(n-l) + (n-l)] 

+ U 2 {N - n) - eo{N - n - 1). 

(4) 

Although it is difficult to obtain exact analytic solutions of Eq. (|3]), we can approximately 
study some interesting phenomena in the high-frequency regime with a; 3> n. To that 
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end, we introduce the function transformation 

C*ni ,n2,n3 (^) '®ni,n2,n3 (^) ^ [ ^ni,n2,n3(t )dt], (5) 

JO 

Then, Eq. ([3]) is transformed to the coupled equations in terms of new amplitudes 
^n\,n2,nzify 


*-BAr,0,0 — ^O-SAf-1,1,0 

X exp[—— F^)t + i— sin(c(;t)] 

(jJ 

]\^ — KyiB ]\;[ — 

X exp[-i(F^+i - F^t + i— sin(a;t)] 

UJ 

“t“ KYi—\B]\j—fi+l^n —\^0 
X exp[-i(F„^_i - F^)t - i— sin(a;t)] 

(jJ 

-|- l^nB]sj—n^n—1,1 

X exp[-i(F,^ - Fl)t + i— sin(ci;t)] 

UJ 

(n=l,2,...,iV-l) 

^-Bo,Af,0 = f^N-lBi^N-1,0 

X exp[—— Fli)t — i— sin(a;t)] 

UJ 

X exp[—— Fl^)t + i— sin(a;t)] 

UJ 

^-®0,n,W—n ^n-®0,n+l,7V—n—1 

X exp[-i(F2^i - Fl)t - i— sin(a;t)] 

UJ 

+ ^n-l-Bo ,n—IjW—n+1 

X exp[-z(F2_^ - Fl)t + i— sin(a;t)] 

UJ 

T ^nBl^n—l,N—n 

X exp[-z(F^ - - i— sin(a;t)] 

UJ 

(n = l,2,...,iV-l) 

*-Bo,0,iV = ^0-5o,l,Ar-l 

X exp[—— Fg^)t — i— sin(a;t)] 


( 6 ) 

To produce the directed motion of dipolar bosons, we assume that the following 
resonance conditions are satished 

(t/o - 17i) = (t/i - U 2 ) = UJ, (7) 

£0 =m = 0, ±1, ±2,... (8) 

Under such resonant conditions, the system will exchange energy with the driving held 
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to bridge the energy gap resulting from strong interactions and static tilt, and the energy 
scale of the system becomes characterized by the natural tunneling coefficient v. When 
the shaking frequency is much larger than the energy scale n, the amount of change in 
^ni,n 2 ,n 3 {t) during a period, T = 27i/u;, can be regarded as being inhnitesimal. Thus Eq. 
(|6]) can be integrated approximately over a period 27i/u by supposing that Bn 3 ,n 2 ,n 3 {t) 
are constants. By averaging the rapidly oscillating exponential terms in Eq. (|6]), we 
obtain a group of approximate equations for the evolution of the amplitudes 


/ \ 

= Koi?Ar-l,l,Oi7-(Ar-l)-m( — ) 

(jJ 


£l. 




iB, 


OJ 

6’l 

“ 1 “ -(N—2n-\-l)—m\ ) 

CJ 

(n = l,2,...,iV-l) 

0,N,0 = I^N-lBi N-l,0>7(N-l)-m{. — ) 

OJ 


Ul 


£l s 


+ I^N-lBo^N-l,l> 7 -(N-l)-m{. — ) 


LJ 
El. 


'i'BQ^n,N—n ^nBQ^ri+l,N—n—l'U{N—2n—l)—m{^ ) 


Ul 

6^1 

“ 1 “ ^n—lBQ rL—l,N—n-\-lJ^(N—2n-\-l)—m\ ) 

CJ 

/ \ 

(n=l,2,...,iV-l) 

• / '^1 \ 

*-Bo,0,V = «^0-Bo,l,Af-li7(Af-l)-m( — ) 

OJ 


Ul 


(9) 


where J7n is the nth-order Bessel function of hrst kind. Equation ([9]) is effective in 
description of the tunneling dynamics of the original system for resonance driven case, 
which is the basis of the following analysis. 


3. Directed tunneling of dipolar bosons under preestablished condition 

Generally, it is hard to solve analytically the large numbers of coupled equations in 
(j9|l. as the dimension of the Hilbert space increases sharply with N. In what follows, 
we are interested in the special case J7k(v-i)-m(^) = 0 or j7(Ar_i)_m(^) = 0, where the 
coupling among the equations in ([9]) is partly removed so that some of the equations 
become closed and analytical. We now proceed to illustrate how to control the tunneling 
processes of a precisely dehned number of bosons. 

When j7k(Ar_i)_m(^) = 0 is selected, the system dynamics is limited in a subspace 
spanned by states \N — 1,1, 0), \N — 2, 2, 0),..., \N — n, n, 0),..., |0, N, 0), and the motion 



7 


of equation ([2]) becomes 


6^1 \ 

iBjsj—Yi ^nB]S[—n—l^n-\-l^oJ^—(N—2n—l)—ra\ ) 

UJ 




l-Sw—n+l,n—1,0*-7^(W—2n+l)—m ( ^ ) 

(n=l,2,...,iV-l) 


iB< 


0,N,0 


6^1 

UJ 


( 10 ) 


It can be learned from Eq. (1T0|1 that the tunneling pathway between wells 2 and 3 is 
shut off and only tunneling between wells 1 and 2 is allowed. 

On the other hand, applying = 0 to Eq. ([9]) yields 




. 




U 


,Sl- 


^uBq n+1 N—n—li^(N—2n—l)—m( ) 

UJ 

6^1 

“ 1 “ /^n—lBQ n—l,N—n+liT(N—2n+l)—m{ ) 

UJ 
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In this case, the tunneling between wells 1 and 2 is prohibited and the tunneling passage 
between wells 2 and 3 is switched on. So far, it has been shown, under the selective 
CDT conditions 77_(Ar_i)_m(^) = 0 or j7(Ar_i)_m(^) = 0, the driven three-site model 
can be reduced to an effective double-well model (ITOH or ffTTj) respectively, in which only 
tunneling between two adjacent wells is allowed. Next, we will present the underlying 
physics behind this decoupling. 

Suppose there is a Fock state \N — i, i, 0), describing that N bosons occupy the left 
and central wells with the right well empty. As this Fock state changes from \N — 
to IA^—i, i—l, 1) (hence one particle is released to the right well), there is a corresponding 
loss of energy 

AEi = [N{U, - U 2 ) - {Uo - U,)] + So. (12) 

Here, the energy loss AEi comes from two parts: one is the loss of interaction energy, 
the other is the energy difference due to the static tilt of triple-well potential. When 
the frequency of driving is chosen such that 


kuj = AEi, (13) 

where k is an integer, the energy of k photons bridges the energy gap between states 
\N — i,i,0) and |A — — 1,1), and, as a result, the tunneling contact disabled by 

both static tilt and strong interaction can be restored in general. This is analogous to 
photon-assisted tunneling. The energy match of Eq. flT^ is realized under the resonance 
conditions (I7])-(|H]). In the fast modulation uj ^ v regime, the effective tunneling 
coefficient between resonant states \N — i,i,0) and \N — i,i — 1,1) is approximately 
renormalized by Jk{ei/u). However, at particular values of the amplitude of driving 


field, Jk{ei/uj) = j7'(Ar_i)+m(^) = = 0, CDT will occur and the tunneling 

dynamics to the right well will be frozen. 

Under such circumstance, the system dynamics is limited in the left and central 
wells, which is governed by the effective motion of equation flTOj) . From Eq. fllOl) we know 
that the effective tunneling coefficient between states |iV — i, i, 0) and |iV — z +1, z — 1, 0) 
is rescaled by a factor of j7L(Ar_2i+i)-m(^)- This renormalization of tunneling coefficient 
can also be interpreted as multiphoton resonances between states |A^ — z,z,0) and 
|A^ — z + l,z — 1,0). As such, if the additional condition 

J-{N-2i+l)-m{ —) = 0 (14) 

OJ 

is satished, no more particle is allowed to tunnel from the central to the left well (hence 
the transition |A^ —z, z, 0) |A^ —z + 1, z —1, 0) becomes prohibited) since the central well 

has already released N — i particles to the left well. Therefore, a desired and prescribed 
number — z of particles can be allowed to tunnel from the central to the left well 
provided that the CDT condition j7L(Ar_i)_m(^) = 0 and the condition of Eq. (ITT]) are 
simultaneously satished. Such two conditions can be achieved by setting 

J,.,{-) = 0,eo = -{N-t)uj. (15) 

LO 

Up to now, we have established the conditions for controlling a dehnite number 
of bosons allowed to tunneling from the central to the left well from an energetics 
argument. The main point of such selective control of tunneling processes is that at 
multiphoton resonances, the system undergoes simultaneously two kinds of CDT effects 
(one is the decoupling of the left-center wells from the right, and the other is particle- 
number-dependent CDT in the left-center wells). 

In the meantime, the decoupling (right-center from left) and the particle-number- 
dependent CDT in the right-center wells can be achieved in a similar manner. The 
particle number dependence of the tunneling coefficient between the central and left 
wells is different from that of the tunneling coefficient between the central and right 
wells. This discrepancy is due to the exactly opposite values of loss of energy resulting 
from constant tilt between the transitions of |A^ — hbO) —)■ |A^ — z -f l,z — 1,0) and 
|0, z, — z) —)■ |0, z — 1, A^ — z -f 1), in spite of the same interaction energy loss. Thus, 
we can switch the directed tunneling of a precisely dehned number of bosons along the 
center-left pathway to the center-right pathway, by only reversing the constant tilt Eq. 
The fact that static tilt breaks the inversion symmetry in the triple well allows for a 
selective control of tunneling processes along different pathways and different directions. 

Following the analysis mentioned above, we summarize the main results as follows: 

(i) The direct tunneling of a dehnite number A^ — z of particles from the central 
well to the left or to the right can be realized under the preestablished conditions 

Ji-i{—) = 0,£o = T{N-i)u. (16) 

OJ 

Here “ — ” is for the directed tunneling along the pathway from the central to the left 
well, and “ -|- ” for the pathway from the central to the right well. 
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(ii) For the initial states |A^ —1,1,0) (|0,1, iV —1)), which may be prepared through 
many-body state engineering using measurements and hxed unitary dynamics [34]. the 
conditions for controlling a desired number i — 1 bosons allowed to tunneling from the 
left (right) to the central well is given by 

J)(^) = 0,£o = T(iV-*-l)w, (17) 

u 

where “ — ” in the second equality corresponds to the direct tunneling from the left to 
the central well, “ + ” to the direct tunneling from the right to the central well. 

(hi) Additionally, if all the Bessel functions in Eq. ffTOj) or flTT]) are non-vanishing, 
N — 1 bosons participate in the tunneling process from the central well to the left 
or right. Thus we can give alternative conditions for such two tunneling processes as 
= 0,eo = muj ^ -{N - i)u, {i ^ 1) and J7(7v_i)-m(^) = ^^6^ = muj ^ 
[N — %)uj^ (i 7^ 1), respectively. 

4. Numerical experiments 

We have checked our theoretical predictions by direct numerical simulations of 
Eq. ([T]). In all simulations, we typically assumed N = 4 bosons and Uq = 75, f/i = 
40, U 2 = 5,u = 35,v = 1, for which the resonance condition ([7]) is satished. 



Figure 1. (color online) Averages = 1,2,3 for different driving parameters: 

(a) Eq = —2a;, Ei/uj = 3.8317; (b) Eq = —co, Eijui = 5.1356; (c) Eq = 0, 
Ei/uj = 6.3802. The initial state is |3,1,0) and the other parameters are chosen as 
Uo = 75, hi = 40 , 1/2 = 5, a; = 35, w = 1. The three cases of Figs. [T](a)-(c) correspond 
to the directed tunneling of zero, one, and two bosons from the left to central well, 
respectively. 
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(a) (b) 



(c) 




Figure 2. (color online) Time evolution of the probability distributions for system 
0 for different driving parameters: (a) eg = = 5.1356; (b) Eq = 2w, 

Ei/uj = 7.0156; (c) eq = — w, Ei/uj = 7.5883. TV = 4 bosons are initially prepared 
in the central well. The other parameters are Uq = 75, Ui = 40, t /2 = 5, u; = 35, w = 1. 
The three cases of Figs.[2Ka)-(c) represent |0,4, 0) —>■ |0, 3,1), |0,4,0) —>■ (70,3,i|0,3,1) + 
(7o,2,2|0,2,2), and |0,4,0) (7o,3,i|0,3,1) + (7o,2,2|0,2,2) + (7o,i,3|0,1,3), respectively. 

As an example, Figs. [I])a)-(c) show the evolution of (rij) = ('0(t)|ajai|'?/l(t)), 
numerically computed from Eq. ([1]) for the initial state |3,1,0) and for three different 
values of driving parameters satisfying conditions (ITTll (therein the second condition 
having “ — ” sign has been applied). In the first case for Eq = —2u [eq = —{N — i — l ) a ;, 
i = 1] and Ei/u = 3.8317 [i7i(^) = Ji{^) = 0], all (uj) maintain their initial values, 
demonstrating that no particle is allowed to tunnel. In the second case for £0 = 

[^0 = —{N — i — 1 ) 0 ;, i = 2] and Ei/u = 5.1356 = J 2 {^) = 0], (ni) oscillates 

between 3.0 and 2.0, and (ns) remains negligible at all times, which demonstrate that 
one particle is allowed to tunnel from the left to right well. Similarly, in the third case 
for £0 = 0 [^0 = —{N — i — l)u, i = 3] and Ei/u = 6.3802 = 0]) we see 

that (ni) oscillates between 3.0 and 1.0, (n 2 ) oscillates between 1.0 and 3.0, and (ns) 
is still zero, which show that two particles are allowed to tunnel from the left to right 
well. These numerical results verify firmly the occurrence of direct tunneling of dipolar 
bosons in a prescribed number under preestablished conditions. 

As demonstrated in Section ([3]), for all particles initially occupying the central well, 
the tunneling of iV — z bosons from the central to right well occurs when the relations 
= 0 and Eq = {N — i)u are satished, and the tunneling oi N — i bosons occurs 
from the central to left well when the condition = 0 still holds and the constant 

tilt £0 is switched to —{N — i)uj. To confirm these predictions, we calculate numerically 
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time evolution of the probability distribution Pni,n 2 ,n 3 = |Cni,n 2 ,n 3 p for the system ([T]) 
with the initial state |0,4, 0). Three sets of our results are shown in Figs.[2]J^a)-(c). In the 
hrst set for Eq = OJ [^o = {N — i)uj, i = 3] and Ei/uj = 5.1356 [j7i-i(^) = = 0], we 

see that full transition between states |0,4, 0) and |0,3,1) occurs, without tunneling to 
other states. This means that only one of N dipolar bosons is allowed to be transferred 
from well 2 to well 3. In the second set for Eq = 2u [^o = {N — i)u, i = 2] and 
Ei/u = 7.0156 [J7i_i(^) = Ji{^) = 0], we see that only transition between |0,4,0) 
and Co,3,i|0, 3,1) + C'o,2,2|0, 2, 2) are allowed, in which two bosons participate in the 
tunneling process along the path between wells 2 and 3. In the third set for eg = 

[eo ^ {N - i)uj, i = 0,1,2] and Ei/uj = 7.5883 = Ji{^) = 0], it can 

be seen that the system experiences transition between state |0,4, 0) and superposition 
state ^ 2^=1 Co,N-n,n\0, N — n,n) with zero population at other states, in which only 
the pathway between wells 2 and 3 is switched on and iV — 1 bosons participate in the 
tunneling process along this tunneling path. 

Moreover, we have demonstrated some other situations by direct numerical 
simulations of Eq. ([1]). In Ref. [33], the authors have numerically exhibited the directed 
tunneling of one particle from the central to right well and directed tunneling of iV — 1 
particles from the central to left well, which can be viewed as two explicit examples in 
our work. Our primary purpose here is to provide a method for precise control of the 
tunneling of a priori prescribed number of dipolar bosons along different pathways and 
along different directions, which will greatly facilitate the control of quantum states. 
The numerical results (not shown) demonstrate that our theoretical predictions are 
still applicable even when the interactions and the driving frequency are not very large 
compared to tunneling rate, indicating that our proposal is more realistic than what it 
seems. 

Before concluding, we present some remarks on our theoretical predictions. 
Dissipation like particle loss presents a major obstacle for long-time coherent control 
of quantum states. The lifetime of bosonic system is principally limited by dissipative 
three-body interactions, which reduces rapidly with the decrease of particle number. 
Like in the schemes of Refs. [25] and [26], high-frequency approximation is valid 
when driving frequency u is of the same order of magnitude or larger than the 
coupling constants Kn,i^n, which requires relatively small numbers of bosons N. 
Presently, moderately small systems of 2-10 atoms can be precisely prepared with high 
purity [35l [36l ET] and long lifetime up to a few tenths of a second[3H], indicating that 
coherent directed tunneling processes should be observable in our considered system. 
Our proposal requires that three interaction parameters (on-site Uq, nearest-neighbor Ui, 
and next-nearest neighbor U 2 ) are subjected to the resonance condition (I7|). According 
to the theoretical calculation in Ref. [5] , the ratio of the nearest-neighbor to next-nearest- 
neighbor interaction U 1 /U 2 depends on the geometry of potential and varies from 4 to 
8 . Note that on-site interaction Uq results from short-range interaction and DDI, in 
which contact interaction between particles can be very precisely controlled by means 
of Feshbach resonance [39], and DDI may be manipulated by varying the shape of a 
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dipolar BEC, the dipole polarization axis, and the trapping geometryWe have 
numerically simulated our main Endings by extending the rigorous resonance condition 
([7]) to nonresonance case Uq — Ui = Ui — U 2 + /3. The numerical results show that our 
theoretical predictions are relatively tolerant against moderate changes in detuning (5. 
As an example, Fig. |3]shows that tunneling process of Fig. |2^a) is insensitive to detuning 
/5. The results imply it is easier to realize experimentally these directed selective- 
tunneling eEects. With the constantly advancing lattice shaking techniques|13]-[22]. 
we expect our results can be tested in the realistic experimental setups. 



Figure 3. (color online) Evolution of the probability distributions versus time for a 
moderate value of detuning /3 = 5. The parameters are the same as those in Fig. [2fa) 
except for next-nearest-neighbor interaction U 2 = 10(d = 5). 


5. Conclusions 

In summary, we have theoretically studied a generalization of many-body selective 
CDT, which enables one to control a priori prescribed number of dipolar bosons 
allowed to tunnel in shaken triple-well potentials. In the high-frequency regimes and 
under the resonance conditions, through rotating-wave (or high-frequency averaging) 
approximation method, we obtain a group of efiectively coupled equations for the 
evolution of probability amplitudes. By adjusting the driving parameters, we can 
decouple these coupled equations and thus establish the conditions for directed tunneling 
of a priori prescribed number of dipolar bosons. Under the preestablished conditions, we 
can transport a desired number of dipolar bosons along difierent pathways and difierent 
directions. 

We expect that our findings can be extended to other multi-well systems and thus 
give us a deep insight into the tunneling dynamics of dipolar condensates in optical 
lattices. Our results may provide an opportunity to manipulate the tunneling of an 
array of dipolar bosons and may be useful for efficient quantum information processing 
and atomic device designing. 
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